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Abstract 

We study the g-state Potts antiferromagnet with g = 3 on the honeycomb lattice. Using an 
analytic argument together with a Monte Carlo simulation, we conclude that this model is 
disordered for all T > 0. We also calculate the ground state entropy to be So/ks = 0.507(10) 
and discuss this result. 
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The effect of ground state disorder and associated nonzero ground state entropy 5*0 has 
been a subject of longstanding interest. A physical example is ice, for which 5*0 = 0.82 ±0.05 
cal/(K-mole), i.e., So/ks = 0.41 ± 0.03 [|l], HQ. Among spin models, an example is the 
Ising antiferromagnet (AF) on the triangular lattice. In the context of this model, Wannier 
argued that a nonzero ground state (g.s.) entropy implies the absence of long-range order, 
viz., staggered magnetization Mst for T > 0. Another example is the Ising AF on the 
kagome lattice @, |^. In both of these Ising models, the nonzero g.s. entropy has the effect 
of removing a phase transition at finite temperature. The Ising AF on the triangular lattice 
is critical at T = 0, while on the kagome lattice, with a larger value of 5*0, it is disordered 
even at T = In these two cases, the nonzero g.s. entropy is associated with frustration. 
However, there are also spin models, such as the antiferromagnetic g-state Potts model 
0-0] on the square (sq) and honeycomb (he) lattice, which exhibit g.s. entropy without 
frustration. Because of the absence of frustration, these models constitute ideally simple 
cases where one can study the effects of ground state entropy on the thermodynamics of a 
statistical mechanical model.0 In contrast to the ferromagnetic (FM) Potts model, which has 
a finite-temperature phase transition for dimensionality d > 1, the question of whether the 
g-state Potts AF has a phase transition at finite (or zero) temperature is more delicate and 
depends on both the value of q and the type of lattice. The q = 3 Potts AF on the square 
lattice has been well studied; an exact result of Baxter showed that it is critical at T = 



111 , in agreement with a renormalization group argument |]T2[, and several Monte Carlo 



simulations have been performed on it |jT3|, |T4|. However, to our knowledge, the behavior of 



the q = 3 Potts AF on the honeycomb lattice has not been definitely established. We report 
here the results of a study of this model. 

The (isotropic, nearest-neighbor, zero-field) g-state Potts model on a lattice A is defined 
by the partition function Z = J2{an} with the Hamiltonian 

^ = E (1) 

(nn') 

where (T„ = 1, g are Z^- valued variables on each site n G A, /? = T~^, and J < for the 
AF case. We define K = pj, a = , x = {a — 1)/^^, and the reduced free energy (per 
site) / = —PF = limTv^oo \iaZ, where denotes the number of sites in the lattice. We 
consider A = he here. 

We first observe that for the g = 2 (Ising) case, the paramagnetic-ferromagnetic (PM-FM) 
^Henceforth, we shall use units such that fcs = 1. 

^Ground state entropy without frustration can also occur in models with continuous variables and inter- 
actions [0. A yet more complicated case is that of quenched disorder with frustration, as in spin glasses. 
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and PM-AFM critical points are both determined by the equation 



15 



./q + 3x-x-^ = (2) 

These are ac = 2 + V3 (PM-FM) and o^af = a"^ = 2-^3 (PM-AFM). (The third root of eq. 
(0) is a complex-temperature singular point a.t a = —1 = ag.) An equivalent representation 
of the partition function is Z = X^g'cg 10) 113' Hi' where G' denotes a subgraph 



of G = A, f = (a — 1), b{G') is the number of bonds and n{G') the number of connected 
components of G' . This enables one to analytically continue the model from positive integral 
q to real q |jl^, ||. Carrying out this analytic continuation and analyzing eq. (|) for the 
critical points, one sees that the points ac{q) and ac,AF('?) increase and decrease, respectively, 
reflecting the fact that as q increases, one must go to lower temperature to achieve FM and 
AFM long-range order. As q reaches the value = (3 + v^)/2 = 2.618.., a^AF decreases 
to 0, i.e., the AFM phase is squeezed out, and there is no longer any finite-temperature 
AF critical point, which now occurs only at T = 0. Note that q^ = = 1 + t, where 
Br = 4cos^(7r/r) is the r'th Beraha number [T^ and r is the golden mean. For q > qz, ac,AF 



is negative, i.e., an unphysical, complex-temperature (CT) singular point. It follows that 
for q > qz and, in particular, for q = 3, the he Potts AF has no critical point or associated 
continuous phase transition at any T > 0.| As q increases from qz to q = 3, the root of 
eq. d^) which, for q < qz, was the PM-AFM critical point ac^ApiQ), nioves leftward from 
the origin. Since for q > qz, there is no longer any physical AFM phase, we shall denote 
this point as ac2', it moves from ac2{qz) = leftward to ac2(3) = —0.1848... Meanwhile, as q 
increases from 2 to 3, (i) the PM-FM critical point adq) moves to the right, through a^qz) = 



(l/2)(3 + \/l5 + 6^5) = 4.1654.. to 0^(3) = 4.4115.., and the disordered, paramagnetic (PM) 
phase (and its complex-temperature generalization) expands accordingly; (ii) the root as{q) 



of (D moves leftward, from -1, through a^iqz) = (l/2)(3 - ^15 + 675) = -1.1654.., to 
a,(3) = -1.2267.. . 

In particular, this argument by analytic continuation in q excludes the possibility, for the 
he lattice, of a massless low-temperature phase with algebraic asymptotic decay of correlation 
functions of the type discussed in Ref. However, this leaves open the possibihty 

that the model might have a first-order transition (with finite correlation length, and hence 

•^From these exact results, recalling the connection Z{q,A,K — —oo) — PA{q), where Pg{<i) is the 
chromatic polynomial for the graph G (lO), we would expect that the behavior of Phc{q) (where A = he 
denotes the thermodynamic limit of the he lattice) would differ for q < and q > q^, a,s would follow if 
the zeros of Phcil) formed a boundary curve in the complex q plane which crossed the real axis at qz and 
separated the regions which include the segments q < qz and q > qz- Given the observation that the crossing 
point of a boundary curve increases by about Aq ~ 0.4 from an 8 x 8 triangular lattice with cylindrical 
boundary conditions (CBC's) to the thermodynamic limit |2^, our expectation is consistent with the finding 
pO| that there is a crossing curve on the 8 x 8 he lattice with CBC's at g ~ 2.2. 
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noncritical). Indeed, this is what did happen for the g = 3 Potts AF on the triangular 
lattice |2^, although for that case, the ground state is only finitely degenerate, so that 
5*0 = 0, and Mgt is nonzero below the transition. In contrast, given that 5*0 is nonzero in the 
present case of the he lattice, the Wannier argument implies that M^t = for all T > 0, so 
that the discontinuity at such a hypothetical transition would have to occur in U but not in 
Mgt- We consider this to render such a transition unlikely but do not know of a proof which 
precludes a discontinuity in short-range order (which enters into U) while long range order, 
Msti remains zero. 

An effective way to study this possibility of a phase transition is to perform Monte Carlo 
simulations of the model, and we have done this.^ For the Monte Carlo simulation, we 
have used two different algorithms to update the spins: the Metropolis algorithm and the 
Swendsen-Wang cluster algorithm (SWCA) |2^. Since the SWCA reduces critical slowing 
down in simulations of models exhibiting a critical points with divergent correlation lengths, 
the agreement of the results obtained from these two algorithms serves as a confirmation of 
our conclusion from the analytic argument above that the model is not critical at T = 0. 
For Metropolis, we used lattices with periodic boundary conditions (BC's) of sizes ranging 
from 8 X 8 to 40 X 40. For SWCA, following Ref. HH, we used lattices with helical BC's 
in the horizontal direction and free BC's in the vertical direction (where the he lattice is 
represented as a brick lattice with horizontal bricks) with sizes ranging from 19 x 19 to 
39 X 39. Typically, we ran for several thousand sweeps through the lattice for thermalization 
before calculating averages. Each average was calculated from 10,000 sweeps through the 
lattice. The full data was obtained as a thermal loop, to test for any hysteresis associated 
with either critical slowing down or metastability. No such hysteresis was observed. The 
results obtained with these two different algorithms were in excellent agreement, differing at 
most by only about 1 %. In Fig. 1 we show measurements of the internal energy per site, 
[/, from the SWCA simulation. 

The intercept and slope at = follow from the high-temperature expansion —U / J = 
{g/2) 1/q + ((g — l)/q^)K + 0{K'^) , where g is the coordination number of the lattice, so 
U/\J\ = 1/2 + {l/3)K + 0{K^) for the g = 3 Potts AF on the he lattice. Clearly, at T = 0, 
i.e., K = —oo, the spins on each bond must be different, so U = —{g/2){6a„an') = 0. As 

"^Two other methods would be (i) to analyze high-temperature series expansions (a first-order transition 
could manifest itself in peculiar behavior of exponents, as in p3|); and (ii) to calculate complex-temperature 
zeros of the partition function and search for a new phase boundary which crosses the positive real a axis 
at a point not included in the roots of (|^). (This relies on the fact that in the thermodynamic limit, these 
zeros typically merge to form curves which separate complex-temperature generalizations of phases. If such 
a boundary were found, the density of zeros near the real axis would yield the critical exponent a, and a = 1 
would suggest a first-order transition.) While these methods are of interest in their own right, we were able 
to obtain convincing evidence for our conclusion from the Monte Carlo method alone. 
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Figure 1: Measurements of internal energy U, as a function oi K — (3J, for the g = 3 Potts 
antiferromagnet on the honeycomb lattice. See text for details. 
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a check on our program, we have also simulated the J > model and obtained excellent 
agreement with the PM-FM phase transition known from eq. (0) to occur at = 1.484. 
Evidently, the data smoothly curves down from the K = value toward the K = —oo value 
as K decreases to —5; there is no indication of any phase transition, in particular, a first- 
order one. The absence of any critical slowing-down for large negative K is in agreement 
with our analytic argument from eq. (0) that the model is not critical at T = 0. The fact 
that the data for the 19 x 19 and 39 x 39 lattices are very close to each other (as was also 
true for the intermediate sizes that we used) shows that it is not necessary to go to larger 
lattice sizes; the present ones are adequate for our conclusion. Indeed, this is not surprising, 
in view of our result that the lattice is disordered for T > (if there had been any indication 
of critical behavior as signalled, e.g. by critical slowing-down, then we would also have run 
simulations on larger lattices). 

Our results imply that the q = 3 Potts AF on the he lattice has the property that in the 
complex a plane, the points a = 1 {K = 0) and a = {K = — oo) are analytically connected 
and a = does not lie on a complex-temperature (CT) phase boundary. We have calculated 
CT zeros of Z on small lattices with periodic boundary conditions and have obtained results 
which are consistent with this conclusion.^ We note that previous studies of the CT zeros of 
the square-lattice Potts model for g = 3 and 4 have shown that the pattern of zeros in the 
Re{a) < region exhibits a significant dependence on the boundary conditions pSf-p? . 



For g > 4, it has been proved that the Potts AF on the he lattice is disordered for all 



T > |2^. This result is quite consistent with our finding, since increasing q increases the 
disorder in the model. 

As part of our study, we have calculated the g.s. entropy So{A,q) = So{hc,3) of the 
model, using the relation 

S{(3) = S{(3 = 0) + (3U{(3) - f U{(3')d(3' (3) 

Jo 

starting the integration at /3 = with S{P = 0) = Ing for the g-state Potts model. We 
found, as in previous work [^, that this provides a very accurate method for calculating 
5*0. For this we used the Metropolis algorithm with periodic BC's for several lattice sizes. 
Since U{K) rapidly approaches its asymptotic value of as -ft^ decreases past about K = —5 
(see Fig. 1), the RHS of (^ rapidly approaches a constant in this region, enabling one to 
obtain the resultant value of S{f3 = oo) for each lattice size. We then performed a fit to 
this data and extrapolated the result to the thermodynamic limit; the results are shown 
in Fig. 2. As a check, we also carried out the analogous calculations for the g = 3 Potts 

^Calculations of CT zeros for this model are also being performed by A. J. Guttmann and I. Jensen. We 
thank these authors for informing us of their work. 
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Figure 2: Measurements of ground state entropy -So, as a function of lattice size, for the 
q — 3 Potts AF on the honeycomb and square lattices. 
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AF on the square lattice. A fit to the finite-size dependence of our sq lattice data agrees 
very well with the form found in Refs. [jl3| and |Q, So{sq, 3) = So{L; sq, 3) + Cgg^s/L'^ with 
Csq,3 = 1.077, and we get So{sq,3) = 0.4317(3), in excellent agreement with the exact value 
So{sq,3) = (3/2)ln(4/3) = 0.4315... [|, |, |]. For the he lattice, as is evident from Fig. 2, 
our measurements do not exhibit the same finite-size dependence as for the sq lattice. An 
empirical function including terms up to yields a good fit to the data (see Fig. 2) and 
gives the L = oo value of the g.s. entropy for the g = 3 he Potts AF 

5o(/ic, 3) = 0.507 ±0.010 (4) 

where the error is an estimate of the uncertainty. This yields W{hc, 3) = 1.66 ± 0.02, where 
W{A,q) = e^°^^''^\ We observe that our results are consistent, to within the uncertainty, 
with the exact expression W{hc, 3) = 5/3. The ratio Rs{-^, q) = 5'(A, q,T = 0)/S{q, T = oo) 
serves as a useful measure of the reduction of disorder in a given model as T decreases from 
oo to 0. Our results yield Rs{hc, 3) = 0.4615 ±0.010 for the g = 3 he Potts AF, which shows 
that the disorder at T = is a substantial fraction of its maximal, T = oo value. 

From the basic relation S = [3U ± / and the property that Ywik^-oo f^U {[3) = 0, as is 
true of the g-state Potts AF models considered here, it follows that 

S'o(A, q) = /(A, q,K = -oo) = lim ln(PA(g)) (5) 

or equivalently W{A, q) = limjv_^oo N^^P\{q). That is, the g.s. entropy is determined by the 
asymptotic behavior of the chromatic polynomial in the thermodynamic limit. Series of the 
form W{A, q) = q({q - l)/g)^^%(A, g), where W{A, g) = 1 ± E^=i WnV"' with y = l/(g - 1) 
were calculated in Ref . |3ll| . It is of interest to compare our result (^ with an estimate from 



the series for W{hc, 3) = 1 + + 2y^^ + 4?/^^ ± calculated through 0{y^^) [0. Because 
of the sign changes in the he series (the coefficients of the first five terms are positive, while 
those of the remaining four terms are negative), it is difficult to make a reliable extrapolation. 
Simply taking the sum yields W{hc, 3) = 1.687, which is agreeably close to our central value, 
1.66. 

In summary, combining analytic arguments and a Monte Carlo simulation, we have 
reached the conclusion that the g = 3 Potts AF on the honeycomb lattice is disordered 
for all T > and have calculated the ground state entropy for this model. 

This research was supported in part by the NSF grant PHY-93-09888. 
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